Bloch Waves and Bloch Bands of Bose-Einstein Condensates in Optical Lattices 



O 

o 

(N 

a: 
c/3 



> 
m 
oo 

O 



c 

o 
o 



13 



Biao Wu, Roberto B. Dicncr, and Qian Niu 
Department of Physics, The University of Texas at Austin, Austin, Texas 78712-1081 

(Dated: February 1, 2008) 

Bloch waves and Bloch band of Bose-Einstein Condensates in optical lattices are studied. We 
provide further evidence for the loop structure in the Bloch band, and compute the critical values 
of the mean-field interaction strength for the Landau instability and the dynamical instability. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db, 73.20.At 



I. INTRODUCTION 

Bose-Einstein condensates (BECs) in optical lattices 
have been attractine, increasing attention from both theo- 
rists [§ §, §, |, 1 and experimentalists §, §. People 
are interested in how the interaction and coherence of this 
system affect the interesting phenomena observed with 
dilute cold atoms in optical lattices , such as Landau- 
Zener tunneling and Bloch oscillations. Recent studies 
have shown that these phenomena are indeed strongly 
influenced by the interaction between atoms. A series of 
novel effects have been discovered, including the nonlin- 
ear Landau-Zener tunneling [Q, the breakdown of Bloch 
oscillations |i| , and dynamical instability |2|, Q . There 
are similar nonlinear periodic systems in other fields, for 
example, the system of the nonlinear guided waves in a 
periodic layered medium |ll| . 

In a one dimensional optical lattice created by two 
counter-propagating off-resonance laser beams, a BEC 
is essentially a one dimensional system when the lateral 
motion can be either neglected H or confined 0. Its 
grand canonical Hamiltonian is 



H=J dx{r(-^+vcosx)^ + -\^-^\ 2 }, 

(1) 

where ip is the macroscopic wave function of the BEC. 
In the above equation, all the variables are scaled to be 
dimensionless with the system's basic parameters, the 
atomic mass to, the wave number ki, of the two laser 
lights, and the average density no of the BEC. The 
strength of the periodic potential v is in units of — — L -, 
the wave function tp in units of y^no, x in units of 



and t in units of ■ 



The coupling constant c = 



4hk 2 L ' —"^""to - fe 2 ) 

where a s > is the s-wave scattering length. A two di- 
mensional version of this system has also received some 
attention rfjjl] . 

In this Brief Report, we study the Bloch bands and 
Bloch waves of a BEC in an optical lattice, and present 
new results that we were unable to obtain in our previous 
studies in Refs.|j], ||]. These new results are possible now 
due primarily to a new development, an exact solution 



found in Ref. Q] . As Bloch bands and Bloch waves are the 
two most important concepts in understanding a linear 
periodic system, they shall also play crucial roles in the 
physics of the nonlinear periodic system (]]]) . Bloch waves 
are the extremum states of the Hamiltonian (|l]) of the 
form 

iP{x,t)=e* k *<p k {x), (2) 

where <j>k{x) is a periodic function of period 2-7r and k 
is the Bloch wave number. Each Bloch wave state (^) 
satisfies the time-independent Gross-Pitaevskii equation 

1 d 

- 2^ + ik) 2 ^ + c\(j}k\ 2 4>k +vcosx<f) k = H<t>k, (3) 

as can be verified by variation of the Hamiltonian (Q). 
Bloch bands are given by the set of eigenenergies [i(k). 




FIG. 1: Schematic drawing of the first and second Bloch 
bands of a BEC in an optical lattice when c > v. 



II. BLOCH BANDS 

In Ref.(I), we studied the tunneling between the two 
lowest bands to see how it is affected by the interac- 
tion. We found that the tunneling is described by a re- 
vised Landau-Zener model, which we call the nonlinear 
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Landau-Zener model. This model predicts a dramatic 
change in the band structure, which is a loop appearing 
at the Brillouin zone edge k = ±1/2 for c/v > 1 (see 
Fig. 0). A direct consequence of this loop structure is the 
breakdown of the Bloch oscillations due to the non-zero 
adiabatic tunneling into the upper band. 

The loop structure is confirmed by an exact solution 
found recently by Bronski et aL(Eq.(10) of Ref.Q), which 
assumes a much simpler form in terms of our notations, 



(4) 



where a± — 
have " 1 



Substituting it into Eq.(pJ), we 



fi — f + c. This solution only exists when c > v, 
and is a Bloch wave at the edge of the Brillouin zone, 
k = 1/2. This Bloch wave carries a non-zero velocity, 

— , while its complex conjugate has an opposite ve- 
This is in sharp contrast with the behavior in 
in which Bloch waves at the 
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locity. 

a linear periodic system, 
zone edge always have zero velocity. This difference con- 
firms the looped band structure. The solution tps and its 
complex conjugate are the two degenerate states at the 
crossing point X (Fig .[[]) . The non-zero velocity carried 
by this Bloch wave is a manifestation of superfluidity of 
BEC. For free particles, the flow e lx ^ 2 is stopped com- 
pletely by Bragg scattering from the periodic potential; 
for the BEC, the flow can no longer be stopped when the 
superfluidity is strong, that is, c > v. 

This loop structure is further supported by our nu- 
merical calculation of the lowest band as shown in 
Fig .||. It is evident that the slope d/z/dfc at the zone edge 
k = ±1/2 becomes non-zero as the interaction strength 
c is increased over the periodic potential strength v, a 
clear indication of the loop structure. However, due to 
the limitation of our numerical method [Q, we are un- 
able to produce directly the loop. An improved numeri- 
cal method is being developed to calculate the loop and 
the higher Bloch bands. 



III. STABILITY OF BLOCH WAVES 

In our second paper we studied the superfluidity 
and stability of the Bloch waves in the lowest band (ex- 
cluding the loop) . We found that the Bloch waves in the 
middle of the Brillouin zone represent super-flows, and 
the other Bloch waves towards the zone edge have both 
a Landau instability and a dynamical instability. More- 
over, we found that these instabilities can disappear from 
all these Bloch waves when the atomic interaction is be- 
yond certain critical values for a fixed lattice strength. 
For easy reference, we call the critical value for the Lan- 
dau instability cl, and the critical value for the dynam- 
ical instability c^. In that work, we were unable to find 
these two critical values because our numerical method 
was not good enough to find accurate Bloch waves at 
the zone edge. Now the exact solution tf^s allows us to 
overcome the difficulty and calculate these two critical 
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FIG. 2: The lowest Bloch bands of of BECs in an optical 
lattice obtained by numerical calculation. Top curve is for 
c = 0.05 > v; middle curve for c = 0.03 = v; bottom curve 
for c = 0.01 < v. The inset is an enlarged version of the tips. 



values, cl and Cd- It is done by studying the stabilities 
of the Bloch wave ipB- Since the Bloch wave at the zone 
edge is the last one to become stable either in terms of 
the Landau instability or dynamically, the critical values 
of c for ipB to become stable are just cl and Cd- 

The physical significance of the two critical values, cl 
and Cd, lies in the way how the Bloch states at k ^ 
are achieved experimentally: the Bloch state at k = is 
first prepared then driven to the desired Bloch states at 
k 7^ by accelerating the optical lattice|l(J. Therefore, 
as the only point connecting the loop to the rest of the 
Bloch band, a stable ipB means that the Bloch states on 
the loop can be accessed and studied experimentally by 
accelerating the optical lattice. 

We first study the Landau instability by analyzing how 
the energy of the system deviates under a small pertur- 
bation. Since the system is periodic, we are allowed to 
write the perturbation as 

yj = ^ B + e** (u(x, q)e l " x + v*(x, q)e- iqx ), (5) 

where q ranges between —1/2 and 1/2, labelling the per- 
turbation mode, and the perturbation functions u and v 
have a periodicity of 2-7T in x. Then the energy deviation 
caused by this perturbation is 



SE 



dx(u*,v*)M{q) 



where 



with 



M(q) = 



C(k) 



_ fC(l/2 + q) 



2 y dx 



ikf 



£(-1/2 + 9) 



v cos x + c 



(6) 



(7) 



(8) 
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and 



■ + Vc s 



2c 



v 

— e 

c 



-2i.r 



2c 



(9) 



If M(q) is positive definite for all —1/2 < q < 1/2, the 
Bloch wave is a local minimum and a super- flow. 
Otherwise, SE can be negative for some q; the Bloch wave 
is a saddle point and has a Landau instability. As already 
noticed in Ref. , the positive definiteness of the matrices 
M(q) for all q's is guaranteed by the positive definiteness 
of M(0). Diagonalizing M(0) for different values of c 
with a fixed v, we obtain the critical value cl, which is 
shown as a dashed line in Fig .||. For the intersection 
point L at v = 0, we have cl = 1/4 = (fc = 1/2) 2 . 




FIG. 3: The critical values of c. The dashed line is cl, the 
critical value of c for all the Bloch waves in the lowest band 
being super-flows; the solid line is Cd, above which all the 
Bloch waves in the lowest band are dynamically stable. 

The dynamical stability of the Bloch wave V's is stud- 
ied by linearizing the Gross-Pitaevskii equation 

dib 1 d 2 ib 
i-jfi = ~Yq^ +cMV + wcosa?V- (10) 

With a procedure similar to the above, we arrive at the 
linearized dynamical equation 



d 

dt \v 



aM(q) 



I 
-I 



(11) 



The dynamical stability is determined by the matrix 
aM(q). If all aM(q) for -1/2 < q < 1/2 have no com- 
plex eigenvalues, then i^b is dynamically stable; other- 
wise, it is unstable. However, as pointed out in Ref.Q, 
the dynamical instability always starts at the perturba- 
tion mode q — 1/2. Therefore, we only need to diagonal- 
ize aM (1/2) to find the critical value Cd- The results are 
shown as the solid line in Fig J|, where the intersection D 
at v = is precisely Cd = 3/16. This lower bound of the 
critical value Cd simply means that when c < 3/16, any 



periodic potential brings the dynamical instability into 
the system. 

The value of Cd at point D is confirmed by analyzing 
the limiting case v -C c, where the matrix aM (1/2) can 
be approximated with a 4 x 4 matrix 



ctM(1/2) 



c-i 

c+ • 

-c 

\ v —c 



\ 
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(12) 



The eigenvalues of this matrix can be found exactly; all of 
them are real only when c > 3/16. Note that the point D 
must be understood in a sense that Cd — >3/16asu— >0 
since precisely at v = the system has no dynamical 
instability. As one may get an impression from Fig.|| 
that the asymptotic behavior of the two curves at large v 
is linear, we want to stress that it is not. Our numerical 
results show that the asymptotic behavior undergoes very 
small oscillations along a straight line, which we have no 
complete understanding. 

Besides the existence of the loop structure shown in 
Fig.|l|, we have not discussed the properties of the states 
on the loop due to the difficulty finding these loop states 
accurately. Here we offer a glimpse of these loop states 
and their properties. Around the crossing point X 
(Fig.0) k = 1/2 + e (|e| < 1), the Bloch wave can be 
approximated to the zeroth oder of e as 



ip(x) 



Vh+l + Vh^l 
2y/h 



ikx 



Vh + l- Vh^l 
2Vh 



i{k—l)x 



(13) 

where h = ^ + ^e. Numerical investigations show that 
most of these loop states are all saddle points which can 
be either dynamically stable or unstable. More numerical 
calculations are needed to further confirm this. 

Finally, we make two remarks. First, the way of defin- 
ing Bloch waves and Bloch bands for the nonlinear sys- 
tem (Q) at the beginning is a natural generalization from 
the linear periodic system. Nevertheless, there is an es- 
sential difference due to the nonlinearity. In the linear 
system (c = 0), the Bloch waves are the only extremum 
states of its Hamiltonian or the only eigenfunctions of 
Eq.(||); for the nonlinear system (|l|), there are possible 
extremum states that are not Bloch waves. 

Second, it is interesting to put the dynamical instabil- 
ity which is discussed in this report and in Ref.|^, [j| 
into perspective. Usually, a quantum dynamics is a reg- 
ular motion because it has discrete eigenvalues thus an 
almost periodic motion no matter its corresponding clas- 
sical dynamics is chaotic or not. In this sense, quantum 
chaos has been called "pseudochaos" jl3j . On contrary, 
the dynamical instability that we have discussed is "true" 
quantum dynamical chaos that deserves more attention 
in the future. 

On the other hand, with the Madelung transformation 
ip(x,t) = p(x,t)e lS ( x,t \ the nonlinear Schrodinger equa- 
tion (fill) can be turned into a set of equations of fluid 
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dynamics. In this regard, the quantum dynamical in- 
stability should be related to the turbulence in the fluid 
dynamics, and we may call it "quantum turbulence" . 



ACKNOWLEDGMENTS 

This work is supported by the NSF, the Robert A. 
Welch Foundation, and the NSF of China. 

REFERENCES 

[1] Biao Wu and Qian Niu, Phys. Rev. A 61 023402 
(2000); O. Zobay and B.M. Garraway, Phys. Rev. A 
61, 033603 (2000) 



[2] Biao Wu and Qian Niu, e-print: |cond-mat /0009455 . 
[3] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86 
2353 (2001) . 

[4] J.C. Bronski, L.D. Carr, B. Deconinck, and J.N. 

Kutz, Phys. Rev. Lett. 86 1402 (2001); 
[5] J.C. Bronski, L.D. Carr, B. Deconinck, J.N. Kutz, 

and K. Promislow, Phys. Rev. E63, 036612 (2001). 



[6] Kirstine Berg-S0rensen and Klaus M0lmer, Phys. 
Rev. A58, 1480(1998); D. Choi and Qian Niu, Phys. 
Rev. Lett. 82 . 

of Opt. B2, 589 (2000). 
and MA. Kasevich, Science 282, 



[9] 



[7] M. Halthaus, J. 
[8] B.P. Anderson 
1686 (1998). 

S. Burger, F.S. Cataliotti, C. Fort, F. Minardi, M. 
Inguscio, M.L. Chiofalo, and M.P. Tosi, Phys. Rev. 
Lett. 86, 4447 (2001). 
[10] C.F. Bharucha, KW. Madison, P.R. Morrow, 
Wilkinson, B. Sundaram, and M.G. Raizen, 
Rev. A 55, R857 (1997); KW. Madison, 
Bharucha, P.R. Morrow, S.R. Wilkinson, Q. Niu, B. 
Sundaram, and M.G. Raizen, App. Phy. B 65, 693 
(1997). 

[11] A. A. Sukhoruk ov and Y.S. Kivshar, e-print: 



S.R. 
Phy. 
C.F. 



nlin.PS/010507S 



[12] B. Deconinck, BA. Frigyik, J.N. Kutz, preprint 

(submitted to Phys. Lett. A.). 
[13] G. Casati, and B. Chirikov, PhysicaD86, 220 (1995). 



